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Chiral phase transitions driven by space-time curvature effects are inves- 

00 . 

, tigated in de Sitter space in the supersymmetric Nambu-Jona-Lasinio model 

q-i with soft supersymmetry breaking. The model is considered to be suitable for 

(2jT), the analysis of possible phase transitions in inflationary universe. It is found 
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Abstract 



that a restoration of the broken chiral symmetry takes place in two patterns 



for increasing curvature : the first order and second order phase transition re- 
spectively depending on initial settings of the four-body interaction parameter 
and the soft supersymmetry breaking parameter. The critical curves express- 
ing the phase boundaries in these parameters are obtained. Cosmological 
implications of the result are discussed in connection with bubble formations 

and the creation of cosmic strings during the inflationary era. 
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In the scenario of the early universe the Higgs mechanism is one of the possible candidates 
in explaining the onset of the inflation era. At the beginning of the inflation era it is assumed 
that the grand unified theory phase is broken down to the quantum chromodynamics and 
electroweak theory phase through the Higgs mechanism. In this connection it is interesting 
to note that the Higgs fields may be composed of some fundamental fermions as in the 
technicolor model and to see the consequence of this idea in the scenario of the inflation. 
On the other hand the supersymmetry is considered to be a vital nature possessed by the 
fundamental unified theory and hence the incorporation of the supersymmetry in composite 
Higgs models seems to be of principal importance. Under these circumstances it is natural 
for us to consider a supersymmetric composite Higgs model in the early stage of the universe 
and to see whether any remarkable effects are drawn during the inflation era. For simplicity 
we adopt the Nambu-Jona-Lasinio (NJL) model as a prototype of the composite Higgs 
model in the present communication. 

At the inflation era the quantum effect of the gravitation is of minor importance while 
the external gravitational field is non- negligible. Hence we are naturally led to the super- 
symmetric NJL model in curved space. Dealing with the composite Higgs fields is essentially 
nonperturbative and does not accept approximate treatments. Accordingly we try to solve 
the problem rigorously working in a specific space-time, the de Sitter space, which possesses 
a maximal symmetry. The de Sitter space is suitable for describing the inflationary universe. 
As a nonperturbative method we rely on the 1/N expansion technique. 

The four-fermion interaction model (which is the basis of the NJL model) in de Sitter 
space has been discussed by several authors and is found to reveal the restoration 

of the broken chiral symmetry for increasing curvature as a second order phase transition. 
The supersymmetric version of the NJL model in curved space was considered by I. L. 
Buchbinder, T. Inagaki and S. D. Odintsov || in the weak curvature limit. They found that 
the chiral symmetry is broken as the curvature increases. Their result is in contrast with 
the result in the nonsupersymmetric NJL model. On the other hand the supersymmetric 
NJL model in the flat space-time has been investigated by several authors [0 in the context 
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of dynamical chiral symmetry breaking. 

We start with the Lagrangian for the supersymmetric Nambu-Jona-Lasinio model in de 
Sitter space expressed in terms of component fields of the superfields 0, 

£ = - V> f V M - p 2 0t _ v> ct V M c - p 2 ct c 

-ip(if - a + i-f 5 7r)ip - ^p 2 , (1) 

where we have kept only terms relevant to the leading order in the 1/N expansion. In Eq. ([I]) 
4> and ip refer to the scalar and spinor component fields of the superfield respectively, N is 
the number of components of these fields, V M the covariant derivative, A the four-fermion 
coupling constant and p the auxiliary scalar field with p 2 = a 2 + ir 2 . We introduce an 
additional term 5£ composed of non-minimal gravitational terms and a soft supersymmetry 
breaking term to the above Lagrangian (JI|) 

5C = -iiR<P ] <P - &R<p c] <p c - A + C V C ), (2) 

where R is the space-time curvature and £i, £2 and A are coupling parameters. We assume 
that the term (||D existed already when the inflation era started. 

The effective potential for the auxiliary field p is calculated in the leading order of the 
1/N expansion such that || 

V(p) = ^— + i[ dsTrS(x,x;s) 
IX Jo 

+2i / dtG(x,x;t) (3) 

J 

where S(x, y; s) and G(x, y; t) represent the fermion and boson propagator in the coordinate 
space with mass s and \ft respectively. The effective potential (|3p is obtained by taking the 
short-distance limit y — > x once the full expressions of these propagators are found in de 
Sitter space. 

The boson propagator in de Sitter space is well-known P-|T2[| and is given for arbitrary 
dimension D by 
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with a± = (-D — 1 ± y (D — l) 2 — Atr 2 ) /2. Here z = a 2 / (2r 2 ) with a the geodesic distance 
and r the radius of de Sitter space which is related to the conventional Hubble parameter 
such that r = 1/H. The scalar propagator for Lagrangian £ + 5C is obtained simply by 
replacing t by t + £-R + A 2 with £ = £x — £2 and R = D(D — l)/r 2 where we have taken 
£1 = £2 for simplicity The fermion propagator is given by 0,0] 

S(x, y; s) = (A(x, y\ s) + B(x, y; s)a.^)U, (5) 

where 

sr 2 ~ D T(a)T(a*) 

A(x,y;s) = 1 



(4tt) d / 2 rp/2 + 1) 

xVT^ 2 F 1 (a,a\D/2 + l;l- z), (6) 

and a = D/2 + isr and {7 is the matrix composed of the Dirac matrices. We do not present 
the explicit expression of the invariant function B(x,y; s) although it is known analytically. 
The reason is that we do not need the explicit expression of this function since the second 
term on the right hand side of Eq. (|5|) disappears for vanishing distance. We find for small 
geodesic distance z ~ that 

TrS(x, x] s) = lim A(x, y; s)TrU 

= Tr[l] lim A(x,y;s), (7) 

with the normalization Till = Tr[l]. Equipped with these propagators for bosons and 
fermions we are now ready to calculate the effective potential @ in an exact form. 

In order to explain our idea in a transparent way we mainly work in 3 space-time dimen- 
sions and give a brief comment on the full extension to 4 dimensions. We use a well-known 
formula which is found in any mathematical table (e.g., WW) to rewrite the boson propagator 
and find for small geodesic distance z ~ 0, 

G(x, y; t) = + -^z/cothvrz/ + 0{z 1 ' 2 ) , (8) 

8nrz 1 ' z Airr 
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where we have defined v = \/tr 2 + 6£ — 1 and ( = £ + A 2 r 2 /6. Note here that we are cal- 
culating the boson propagator for Lagrangian £ + 5C The fermion propagator at vanishing 
distance in 3 dimensions is obtained in a similar way. Using the relation (^) we have 



TrS(x, x; s) = Tr[l] 



is is / 1 2 2 \ tanhzrsr 

lixrz 1 / 2 Aixr \4 / sr 

+0(z 1 / 2 ) 



(9) 



Substituting propagators (H) and (01) into effective potential (BT) we find 



\ tanliTrsr 



sr 



dtvcothnv, (10) 



2nr Jo 

where we have set Tr[l] = 4 (We adopt the reducible representation of the Clifford algebra of 
Dirac matrices to afford the existence of 75). It is important to note here that the divergence 
present in both expressions (^) and cancel out in the effective potential (PP. The origin 
of this cancellation may be traced back to the supersymmetry of our model. 
The gap equation is given by 



— \J p 2 r 2 + Q( — lcoth(7r^/ p 2 r 2 + 6( — 1 

where by V'(p 2 ) we mean the differentiation of the effective potential with respect to p 2 . 
Just by observing the left hand side of the gap equation (fill) we find that it is a function only 
of pr with parameters nr/\ and (. Thus the solution pr of this gap equation is completely 
specified by two parameters nr/X and ( respectively. This fact suggests that the phase 
diagram of the model will be given on a plane specified by these two parameters. 

The direct observation of Eq. (JTTj) shows that it has at most two solutions for pr and so 
the shape of the effective potential (0) is of three types: the symmetric type S and type I 
and II which are shown in Fig. 1. The phase diagram of the model is given by the numerical 
analysis of the effective potential fllCf ) and the gap equation (0), and is given in Fig. 2. 



2nr 



nr (\ 9 r,\ tanli7rpr 

— + - + A 2 ] 



A \4 J pr 

~ =0, (11) 
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The region above the dashed line and the right half of the solid line in Fig. 2 represents a 
broken phase with the effective potential of the shape II given in Fig. 1. The small region 
between the solid line and the dashed line in Fig. 2 corresponds to a broken phase with the 
effective potential of the shape I in Fig. 1. The region below the whole solid line represents 
a symmetric phase where the shape of the effective potential is of the single- well type S. 

The boundaries which separate the above three phases are determined as follows: The 
dashed line and the right half of the solid line constitute the boundary of the region charac- 
terized by the potential of the shape II. The boundary is determined by the condition 

V'{0) = 0. (12) 

The above equation reduces to 



irr 7r 



x ( /6C-lcoth7rV6C-l = 0. (13) 

The condition for the line separating the phases of type I and type S is given by V(p*) = 
where p* = max{p 1; p 2 } with pi and p 2 two solutions of the gap equation (|TT|). 

The branching point C in Fig. 2 is of special interest. It is a critical point which divides 
the broken phase into type I and type II. At the branching point C the following conditions 
are found to be satisfied simultaneously: 

V'(0) = 0, V"(0) = 0. (14) 

The condition (|14D is explicitly given by Eq. fll~3|) and 

1 ^ , 1 _ coth7ry/6C-l 

12 2 sinh 2 vt^C - 1 2ttV6C - 1 ' 

From Eq. fllSD we find that £ = 0.290138(= £*) at the branching point C. By substituting 



this value for £ in Eq. (|T3| ) we obtain nr/X = 0.083063(= 77*) at the branching point C. 

Now let us discuss the time-evolution of the chiral structure of the model assuming that 
the inflationary era is well described by the effective potential in de Sitter space. It is natural 
to assume that the curvature slowly decreases (r increases) as the universe evolves. First 
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let us consider the case where the soft supersymmetry breaking term is not included, i.e., 
A = 0. In this case it is easily seen in Fig. 2 that we move from left to right by increasing 
radius r with ((= £) and A fixed. By direct observation of the gap equation one can easily 
show that, if the parameter £ is kept below £ = 1/4 which is the value of £ corresponding 
to irr/\ = 0, the effective potential stays in the type S as r increases and hence the chiral 
symmetry is preserved (there is no phase transition). If the parameter £ is kept in the region 
1/4 < £ < £*, the effective potential changes its shape from type I to type S as r increases. 
Hence the broken chiral symmetry is restored as r increases through the first order phase 
transition. If the parameter £ is kept above (*, the effective potential changes its shape from 
type II to type S and so the transition is of the second order. Thus in the case with A = the 
chiral symmetry restoration occurs as the universe evolves. This situation reflects the fact 
that the non-minimal gravitational coupling terms, which break the supersymmetry that 
protects the chiral symmetry for any value of A and r, are effective only at large curvature. 
It is also interesting to examine the pattern of the phase transitions when the parameter 
£ is changed with r and A fixed. In this case we easily conclude that the chiral symmetry 
is broken by increasing £ through the first order phase transition for nr/\ < rj* while the 
transition is of the second order for nr/\ > rj*. 

We next focus our attention on the role of the soft supersymmetry breaking terms. Here 
we assume the conformal gravitational coupling in 3 dimension, i.e., £ = 1/8, for simplicity. 
This assumption is not essential for discussions in the following as long as £ < 1/4. In 
Fig. 3 we show trajcetories on the phase diagram as the radius r increases with A and A 
fixed (See dot-dashed lines (a)-(c) in Fig. 3). The trajectories are specified by the equation 
C = 1/8 + (AA/7r) 2 r]/6 with rj = irr/\ for each value of AX/n. Increasing the radius r 
along the curve (a) one experiences a first order phase transition and along the curve (c) 
a second order one. The critical case is shown by the curve (b). These curves show that 
the chiral symmetry breaking can occur in the model with the soft supersymmetry breaking 
terms when the radius r increases or when the curvature of the universe decreases during 
the inflation. By observing the behaviors of these curves we find that the types of the phase 
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transitions are classified as follows: If AX/tt > a(= 29.5yC* — 1/8 = 12.0) is satisfied, the 
phase transition is of the first order. The case 1 < AX/tt < a corresponds to the second 
order phase transition, and the phase transition does not occur if AX/tt < 1 is satisfied. 

In summary we have investigated the supersymmetric NJL model with the non-minimal 
gravitational coupling terms and soft supersymmetry breaking terms in de Sitter space. The 
phase structure is completely clarified in the model of 3 dimensions. We have found that both 
the non-minimal gravitational terms and the soft supersymmetry breaking terms lead to the 
phase transition phenomena as radius r increases or the curvature of the universe decreases. 
These two terms work in different ways. The soft supersymmetry breaking terms lead to the 
chiral symmetry breaking as the curvature decreases while the non-minimal gravitational 
terms lead to the restoration of the chiral symmetry. The extension of our work to the 4 
dimensional case is straightforward although we have to rely on the full numerical estimates. 
The result of the analysis will be presented elsewhere. 

The cosmological consequences of the phase transitions found here seems to be quite in- 
teresting. A production of cosmic strings is expected because the [/(l)-symmetry is broken 
in the model. The formation of topological defects in the inflationary universe has been 
studied in many references (e.g., It is well-known that the cosmic strings are con- 

strained observationally as Gji ~ (9(1CT 5 ) where G is the gravitational constant and \i is 
the line density of the cosmic string [0 . This condition sets a constraint on our model and 
it is an interesting problem to investigate whether the condition is satisfied in the model. 
In order to apply our model to a more realistic inflationary scenario we need to deal with 
the model in 4 dimensions. At the moment we have found that the phase structure in 4 
dimensions is very similar to that in 3 dimensions. The expression of the effective potential 
in 4 dimensions is rather cumbersome and the detailed analysis is still in progress. 

Another comment is directed to the possibility of creation of an open universe in our 
model. Inflation models for an open universe have been investigated by using a bounce 



solution for a bubble formation | 18| . The trajectory (a) in Fig. 3 indicates that bubble 
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nucleations may occur in the region I. It is worthy to investigate whether our model can 
provide a successful model for an open universe or not. 

The authors would like to thank Roberto Camporesi, Atsushi Higuchi, Tomohiro Ina- 
gaki and Misao Sasaki for enlightening discussions and useful correspondences. Two of the 
authors (T. M. and K. Y.) are indebted to Monbusho Fund (Grant-in-Aid for Scientific 
Research (C) from the Ministry of Education, Science and Culture with contract numbers 
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FIGURES 
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FIG. 1. Typical behaviors I and II of the effective potential. 
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FIG. 2. The phase diagram on the parameter plane in Q and irr/X. The point denoted by C is 
the branching point where (tvt/X, C) = (r/*,C*) with rj* = 0.083063 and C* = 0.290138. 
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FIG. 3. Typical trajectories for varying r with A and A fixed, (a) A\/tt = 22.3; 
(b)AA/vr = 12.0; (c) AA/vr = 3.2. Here £ is fixed so that f = 1/8. 
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